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A dynamical model for the dusty ring in the Coalsack 
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Abstract. Lada et al. recently presented a detailed near-infrared extinction map of Globule G2 in the Coalsack molecular cloud 
complex, showing that this starless core has a well-defined central extinction minimum. We propose a model for G2 in which a 
rapid increase in external pressure is driving an approximately symmetric compression wave into the core. The rapid increase 
in external pressure could arise because the core has recently been assimilated by the Coalsack cloud complex, or because the 
Coalsack has recently been created by two large-scale converging flows. The resulting compression wave has not yet converged 
on the centre of the core, so there is a central rarefaction. The compression wave has increased the density in the swept-up gas 
by about a factor of ten, and accelerated it inwards to speeds of order 0.4kms _1 . It is shown that even small levels of initial 
turbulence destroy the ring seen in projection almost completely. In the scenario of strong external compression that we are 
proposing this implies that the initial turbulent energy in this globule is such that E turh /E grav < 2%. Protostar formation should 
occur in about 40, 000 years. 
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1. Introduction 

In the current paradigm for star formation, single stars, bi- 
nary and multiple stars, and even small-N clusters, are pre- 
sumed to form following the collapse and fragmentation of a 
dense molecular core (e.g. Andre, Ward-Thompson & Barsony, 
2000). However, many details of this paradigm remain to be un- 
derstood, in particular the processes which form dense molec- 
ular cores in the first place, and thereby determine the initial 
conditions for collapse and fragmentation. 

In an attempt to understand these processes, Lada and col- 
leagues have explored the structures of several cores, using the 
NICE technique (Alves et al. 1998). These include B68 (Alves, 
Lada & Lada (2001), where there is evidence that the core is 
oscillating, in accordance with the predictions of Matsumoto 
& Hanawa (2003) and Hennebelle (2003), and B335 (Harvey 
et al. 2001), where the density profile suggests that the core 
has already become unstable against collapse. In a more re- 
cent application of the NICE technique, Lada et al. (2004) have 
mapped Tapia's Globule 2 in the Coalsack (Tapia, 1973). 

The Coalsack is an extended cloud complex at a distance 
D ~ 150 to 175 pc (Cambresy, 1999; Rodgers 1960). If we 
adopt the smaller distance, the angular extent of the complex 
(~ 6°) corresponds to a linear extent of ~ 15 pc, and the mass of 
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molecular gas is ~ 3, 000 M (Nyman, Bronfman & Thaddeus, 
1989). The highly structured molecular gas mapped by Nyman 
et al. is probably accompanied by an envelope and/or infill of 
somewhat more diffuse HI, providing shielding from the ambi- 
ent H2 -dissociating UV radiation field. This HI would be very 
hard to detect unambiguously, since the systemic velocity of 

the Coalsack is very low ( 5 to kms -1 ) and there would 

therefore be confusion with local 21 cm emitting gas. 

Despite its molecular content, and therefore presumably 
relatively high density, the Coalsack cloud complex shows 
no unambiguous evidence for ongoing star formation. There 
are no flare stars (Weaver 1973); no T Tauri stars (Weaver 
1974a,b); no Ha emission stars (Schwartz 1977); and no IRAS 
sources associated with either the extinction peaks or CO emis- 
sion peaks (Nyman, Bronfman & Thaddeus, 1989). Reipurth 
detected a single HH object, R10, in the direction of the 
Coalsack, and Eaton et al. (1990) showed it to be a bipolar neb- 
ula with a heavily obscured illuminating star. However, Kato et 
al. (1999) question whether R10 is really associated with the 
Coalsack, and it is certainly a long way from Globule 2 (here- 
after G2). 

It is therefore possible that G2, which appears to be the 
densest core in the Coalsack, is in an early stage of condensing 
to form a protostar (or protostars) and, if so, its detailed struc- 
ture could contain important clues about how cores form. To 
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this end, Lada et al. (2004) have obtained deep near-IR pho- 
tometry of ~ 24, 000 stars seen through G2, and have used this 
data to construct an extinction map of the core with unprece- 
dented resolution. 

The most remarkable feature of this map is that the core 
has a local extinction minimum in the middle, and the max- 
imum extinction is in a ring-like structure around this mini- 
mum. Lada et al. discuss some possible explanations for this 
structure, under the general headings: 'flattened geometry and 
magnetic fields' and 'spherical geometry: a core in transition'. 
We here present a detailed model under the second heading, in 
which G2 is a core responding to a sudden increase in external 
pressure as investigated in Hennebelle et al. (2003, 2004). This 
increase in external pressure is driving a compression wave into 
the cloud, producing a dense shell which - seen in projection 
- is visible as a ring of enhanced extinction. The broad 12 C 18 
(J = 2 — > 1) emission line observed by Lada et al. along a 
single line of sight through the core is then attributable to the 
extra trans-sonic velocity generated by the inward-propagating 
compression wave. 

In Section 2 we describe briefly the code used to simulate 
this model, and the initial conditions adopted. In Section 3 we 
present the results obtained and compare them to the obser- 
vations; it should be emphasised at the outset that we do not 
perform an extensive parameter search to obtain the best fit to 
the observational data, since this does not, at the present time, 
seem justified by the scope of the observations. In Section 4 we 
discuss the results, attempt to justify key aspects of the model, 
and summarise our main conclusions. 



2. Numerical method and initial conditions 

The simulations were performed using the SPH code 
DRAGON, with ~ 500, 000 particles. This code invokes an 
oct-tree to find neighbours and calculate gravitational accel- 
erations. The particle smoothing lengths are adjusted so that 
each particle has N NElB =50 + 5 neighbours. A second-order 
Runge-Kutta integration scheme is used, with multiple parti- 
cle time-steps. Other details of the code are given in Goodwin, 
Whitworth & Ward-Thompson (2004). 

In the first instance we present a simulation with zero tur- 
bulence, since this seems to capture well the azimuthally av- 
eraged features of the ring-like structure reported by Lada 
et al. (2005). The initial conditions for this simulation are a 
truncated Bonnor-Ebert sphere with boundary at radius R B = 
3a /(47rGp c ) 1/2 , where a is the isothermal sound speed and p c 
is the central density. We set a = 0.2 kms" 1 (corresponding to 
molecular gas at T Q = 10 K) and p c = 0.2 x 10" 20 g cm -3 (cor- 
responding to n c = 0.5 x 10 3 H2 cm -3 ). This gives a core radius 
R Q = 0.4 pc and a total core mass M TOTAL = 6.5 M Q . The cloud 
is in stable equilibrium with external pressure P EXT = P = 
0.22 x 10 -12 erg cm" 3 = 0.16 x 10 4 k B cm" 3 K. The initial con- 
ditions are set up with a random but settled (i.e. non-crystalline) 
distribution of particles in detailed hydrostatic equilibrium, and 
the subsequent evolution is isothermal (i.e. P = a*p, with uni- 
form and constant a ). At t - 0, we increase P EXT at a constant 



Table 1. Parameters of the initial turbulence in the different 
simulations. 
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rate, by a factor of 10, over a thirtieth of a crossing time, i.e. 

K. t<0; 
P EXT = (1 + 270f/r sc ) P , 0<t< f sc /30 ; (1) 
[l0P , f>f sc /30. 

The sound-crossing time is given by t sc = R B /a = 2Myr, 
and so the new increased external pressure is reached after ~ 
0.07 Myr. 

We then explore how the evolution of the core is modified 
if it has a modest level of internal turbulence. This is moti- 
vated by the fact that the ring-like structure observed by Lada 
et al. presents significant departures from an ideal axisymmet- 
ric ring. We therefore want to explore wheter turbulence is able 
to produce similar patern and for which level of turbulence. 

The turbulence is characterised by two parameters: (i) the 
ratio of turbulent energy, £ XURB , to self-gravitational potential 
energy, D, i.e. Qf TURB = E TUgB /\Q\\ and (ii) the exponent n in the 
turbulent power spectrum, P(k) oc k~ n . The turbulent parame- 
ters in the different simulations are given in Table 1 

3. Simulation results 

3.1. Spherical case 

Fig. 1 shows the profiles of volume-density, n(r), inflow ve- 
locity, w(r), and column-density, N(b), generated by the non- 
turbulent simulation (Simulation 1 in Table |3 at three times, 
t = 1.170, 1.183, and 1.198 Myr (i.e. well after the external 
pressure has reached its new increased value, 10 P Q ). Here r 
is distance from the centre of the core and b is the impact pa- 
rameter of the line of sight relative to the centre of the core. 
The density in the swept-up shell has been increased from 
n ~ 0.3 x 10 3 cm -3 to n ~ 0.6 x 10 4 cm" 3 , and the swept-up 
shell is travelling inwards at u ~ 0.4 kms -1 (i.e. Mach-2). 

If we adopt the relation N = 0.2 x 10 22 cm" 2 x A v (where 
N is the column-density of hydrogen nuclei in all forms, and 
A v is the visual extinction in magnitudes), we can convert the 
azimuthally averaged extinction profile obtained by Lada et al. 
(2004) from their near-IR measurements through G2 (their Fig. 
3) into a column-density profile. The result of this conversion 
is shown as a dashed line in the righthand panels of Fig. 1. 
We see that the observed column-density profile is fitted best 
by the simulated column-density profile at t = 1.183 Myr. The 
simulated profile has a minimum of N ~ 0.16 x 10 23 cm -2 in 
the centre, and a maximum of N ~ 0.24 x 10 23 cm -2 at radii 
r ~ 0.04 pc, close to the maximum of the observed profile. We 
note that, with a spherically symmetric model, it requires a very 
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Fig. 1. The lefthand column shows volume-density profiles, n(r) (where n is in H2 cm -3 and r is in pc. The middle column sows 
radial velocity profiles, u(r) (where u is in kms -1 ) and r is in pc). The righthand column shows column-density profiles Nib) 
(where N is is H2 cm -2 and b is in pc. The simulation results are represented by solid lines, and the column-density profile derived 
from the observed extinction is represented by a dashed line. The top row is at time t = 1.170 Myr; the middle row is at time 
t = 1.183 Myr; and the bottom row is at time t = 1.198 Myr. 



deep - approximately tenfold - drop in the volume density n at 
the centre to produce a ~ 30% dip in the column-density N m 9 
it is very unlikely that the subsonic velocities involved in an 
oscillating core could reproduce this feature. 

It is worth to stress that our model predicts a column density 
in the external part of the cloud (x > 0.1 pc) which is roughly 
two times lower than the observed column-density. This dis- 
crepency could be due to the fact that we do not model the gas 
outside the core which may have, either a direct contribution to 
the total column-density or could have been accreted in part by 
the core therefore enhancing the core density in its outer part. 

Lada et al. (2004) also report a single line profile for the 
12 C 18 (/ = 2 — > 1) transition from a line of sight which is 
displaced from the notional centre of symmetry of the core by 
AS = 68.5" and Aa = 69.", giving a net angular displacement 
of ~ 100". At a distance of D = 150 pc, this corresponds to an 
impact parameter b ~ 0.08 pc. We can therefore predict the line 
profile expected from the simulation, using 
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Fig. 2. 12 C 18 (J = 2 -> 1) line profiles at time t = 1.183 Myr, 
from lines of sight at different impact parameters b, relative 
to the centre of the simulated core. The impact parameters are 
b = 0, 0.04, 0.08 and 0. 12 pc. 



where v is the radial velocity (i.e. frequency) observed, 
j(n(r),T ) is the integrated volume emissivity in the line, s 
is distance along the line of sight (measured from the tangent 

point), r(s) = (b 2 + s 2 ^ 2 , v'(s) = u (r(s)) s/r(s) is the inflow 
velocity projected along the line of sight, and 0(v'-v) is the pro- 



file function (reflecting thermal, turbulent and natural broaden- 
ing of the line). 

Fig. 2 shows emission line profiles for 12 C 18 (/ = 
2 — > 1) on lines of sight with impact parameters b = 
0, 0.04, 0.08 and 0. 12 pc. The simulated line profile at b = 
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Fig. 3. Column density maps from simulations with initial tur- Fig. 4. Same as Fig. |3 for a Tl]RB = 0.02, n = 5/3. Top panel 
bulence. Parameters are: a Tl]RB = 0.01, n = 5/3. Top panel is for is for time equal to 1.17 Myr whereas second panel is for 
time equal to 1.17 Myr whereas second panel is for 1.18 Myr. 1.185 Myr. 



0.08 pc matches the width of the broad line component ob- 
served by Lada et al. (2004) quite well. We concur with their 
interpretation of the narrow line component as due to a small- 
scale turbulent fluctuation, and we would not expect the present 
model to reproduce this. From the profiles presented in Fig. 2, 
it should be straightforward to test the model proposed here by 
measuring 12 C 18 (/ = 2 —> 1) line profiles at different posi- 
tions within the G2 core. Positions within ~ 50" of the centre 
should present double-peaked line profiles, due to the receding 
and approaching sides of the compression wave. 



3.2. Effect of turbulence 

Here we present the results of three simulations in which the 
initial conditions are the same as in Simulation 1 , except that 
a turbulent velocity field is superimposed at the outset. The 
figures all show column density maps. The different turbu- 
lence parameter combinations are given in Table 1. We note 
that the levels of turbulence invoked are very low £# TURB — 
0.01 and 0.02. The turbulent power spectra have exponents 
n = 5/3 (corresponding to Kolmogorov) and (i.e. more power 
on small scales). 

Fig. 3 displays the column densities at two times steps close 
to t= 1.18 Myr which matches best the observation. It shows 
that even with a TURB = 0.01 and n = 5/3 (Simulation 2) the ring 
is significantly distorted. It is still roundish but is fragmented 
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Fig. 5. Same as Fig.[3]for a Tl]RB = 0.02, n = 0. Top panel is for 
time equal to 1.175 Myr whereas second panel is for 1.20 Myr. 



in denser clumps. We note that this is actually the case for the 
ring observed by Lada et al. 

Fig. 4 shows that with Qf TURB = 0.02 and n = 5/3 
(Simulation 3) the ring is lost rapidly and completely. The gas 
rapidly concentrates in one single dense clump. 

The only way to retrieve the ring with a Tl]RB = 0.02 is to 
concentrate turbulent energy at short wavelengths by decreas- 
ing n. This is illustrated in Fig. 5, where Qf TURB =0.02 and n = 
(Simulation 4), and the ring stays almost axisymmetric during 
a longer amount of time. However, such a flat power spectrum 
is unrealistic since usually Kolmogorov-like powerspectra are 
found in observations and in numerical simulations. We note 
that for large values of Of TURB (> 5%), the ring is destroyed even 
with n - 0. 



We conclude that, within the scope of our model, the ring- 
like structure observed by Lada et al. (2005) requires the initial 
core to have a very low level of turbulence. This suggests that 
the initial core has been very settled before it was subjected to 
increased pressure so that there has been enough time for the 
turbulence to decay almost completely. 

4. Discussion and Conclusions 

We have modelled Tapia's Globule G2 in the Coalsack as a core 
which is responding to a rapid increase in external pressure. 
The resulting compression wave has not yet reached the cen- 
tre, and therefore there is a central rarefaction which manifests 
itself as a dip in the extinction towards the centre, as observed 
by Lada et al. (2004). This model requires that the core was al- 
ready in existence as a coherent entity with relatively low mean 
density, n mT ~ 0.4 x 10 3 cm" 3 , and large extent, r INIT ~ 0.4 pc, 
and that the external pressure acting on it increased by a factor 
~ 10 (from ~ 0.16 x 10 4 k B cm" 3 K to ~ 0.16 x 10 5 k B cm" 3 K) 
in a small fraction of a sound-crossing time, t sc /30 ~ 0.07 Myr. 
The advantages of such external triggers for synchronising the 
formation of stars in clusters and generating high multiplicity 
have been emphasised by Whitworth et al. (1996). 

One possibility is that the Coalsack has been involved in, 
or even created by, a collision between two large-scale turbu- 
lent flows. To deliver the required large pressure increase in 
the required short time requires a flow with minimum den- 
sity n MlN ~ 2 cm" 3 and minimum speed v MIN ~ 6kms _1 . 
Neither requirement is extreme. Indeed the models of Audit & 
Hennebelle (2005) show that pressure increases of the magni- 
tude invoked by our model are very common in the atomic gas 
from which molecular complexes like the Coalsack are pre- 
sumed to form. 

A second possibility is that the core has fallen into a denser 
environment and been compressed by the increased ambient 
pressure. The local cloud within which G2 resides is Region 
II in the scheme of Nyman et al. (1989); see also Kato et al., 
(1999). This cloud has mass M ~ 400 M o and cross-sectional 
area A ~ 9pc 2 , so its mean surface density is £ ~ M/A ~ 
0.01 gem -2 , and therefore, assuming that it is self-gravitating, 
its central pressure must be P c ~ G£ 2 ~ 0.5 x 10 5 k B cm -3 K. It 
follows that the core only needs to travel about one third of the 
cloud radius, i.e. ~ 0.5 pc, to experience the required increase 
in external pressure, and it can do this in t sc /30 ~ 0.07 Myr 
by travelling at ~ 7 kms -1 , which is a typical bulk velocity 
for a cool cloud in the local interstellar medium. Since the 
measured radial velocity of G2 is 6 kms -1 , and the mea- 
sured radial velocity of the Coalsack is ~ -5 to kms -1 
(Nyman et al. 1989), the relative radial velocity is in the range 
~ 1 to 5 kms -1 . Thus if the net velocity of G2 relative to 
the Coalsack is ~ 7 kms -1 , the chance of its radial compo- 
nent being in this range is ~ 71%. Given that we have only 
one case, we cannot take statistical arguments any further, but 
this at least illustrates that the model is not based on a highly 
unlikely coincidence. We note that this estimation is based on 
the assumption that the main effect of the globule falling into 
a denser region would be the pressure increasement. However 
other processes like accretion of more material into the original 
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globule could certainly have a significant impact and promote 
the formation of a core out of equilibrium. 

We have also shown that within the scope of our model the 
initial level of turbulence in the core (before t = 0) must have 
been very low (<*turb < 5%) in order for the ring-like structure 
to develop in a coherent manner. Turbulence promotes the frag- 
mentation of this ringlike structure seen in projection and gen- 
erates dense clumps which look like some of the substructures 
observed by Lada et al. in their ringlike structures. However we 
note that unlike in our simulations, in the observations these 
overdensities seem to protrude from the ring. We do not know 
at this stage the reason of this. This may be due to the setting up 
of turbulent fluctuations with random phases or to the incom- 
plete treatment of the physics of the core since here we ignore 
magnetic field and large scale anisotropy. Our numerical res- 
olution could also be insufficient to describe accurately such 
small scale structures. 

We reiterate that we have not performed an extensive pa- 
rameter search to obtain a best-fit model. The observational 
constraints (the azimuthally averaged, radial extinction pro- 
file; and a single, off-centre, approximately pencil-beam, line- 
profile) do not justify detailed modelling, so we are here only 
attempting a feasibility study. The next step in testing the model 
should be to perform further line observations to test the kine- 
matic predictions of the model (see Fig. 2). We hope to obtain 
these observations soon. Another line of investigation which 
should be pursued is to explore numerically the departures from 
spherical symmetry which can be expected because the pres- 
sure increase acting on the core is likely to be anisotropic. 

In the meantime, it is appropriate to contrast the predictions 
of our model with those of the magnetic model developed by 
Li & Nakamura (2002). 

(i) The magnetic model requires that the core be viewed 
from close to the axis of symmetry. In contrast, our model can 
be viewed from any angle, since it is spherically symmetric. 

(ii) The magnetic model requires the magnetic field to be 
mainly along the line of sight. Lada et al. (2004) show that the 
field component in the plane of the sky (estimated from the 
dispersion in polarisation angles of background stars, using the 
Chandraskhar-Fermi equation) is rather close to the equiparti- 
tion value, and they point out that this rather large field com- 
ponent in the plane of the skyis hard to reconcile with the idea 
that the main component is along the line of sight. 

(iii) G2 is located towards the edge of the cloud complex, 
so it is difficult to see how it can evolve quasistatically in the 
ordered way described by Li & Nakamura (2002). This pre- 
sumably requires the field to be anchored in some larger- scale 
structure, but no such structure is evident on the CO maps of 
Nyman et al. (1989) and Kato et al. (1999). 

The model proposed here falls naturally within the cate- 
gory 'spherical geometry: a core in transition' defined by Lada 
e al. (2004). The question then arises as to whether the core 
will eventually settle into an equilibrium state (a Bonnor-Ebert 
sphere) or whether it will continue contracting to form a star 
(or stars) at the centre. Indeed as in Hennebelle et al. (2003, 
2004) later times of the simulations show that when the com- 
pression wave converges on the centre of the core it will create 



a primary protostar. The compression wave is expected to reach 
the centre in about 40, 000 years. 
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